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Abstract 

Classical dynamics in SU(2) Matrix theory is investigated. A classical chaos- 
order transition is found. For the angular momentum small enough (even for small 
coupling constant) the system exhibits a chaotic behavior, for angular momentum 
large enough the system is regular. 
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1 Introduction 



Matrix theory []J is a surprisingly simple quantum mechanical model that is able to 
describe some major properties of superstring theory. Therefore the model obviously 
deserves a thorough study. Calculation of physical quantities is reduced to the appropriate 
calculations in the matrix quantum mechanics. A system of N Dirichlet zero-branes is 
described in terms of nine N x N Hermitian matrices Xi,i = 1, ...,9 together with their 
fermionic superpartners. The action can be regarded as the dimensional reduction of 
ten-dimensional SU(N) supersymmetric Yang-Mills to (0 + 1) space-time dimensions: 

S = J dtTr(^D t X t D t X t + ^[X i: Xj][X u X^) + (fermtons), (1) 
where D t = d t + iA . 

The action ([!]) was considered in the theory of eleven-dimensional supermembranes 
in |5|, [3], |J and in the dynamics of D-particles in |5|, [7J. In the original formulation 
of the conjectured correspondence between M-theory and M(atrix) theory the large N 
limit was assumed. A more recent formulation [§] deals with finite N. In the last year 
the Matrix theory was a subject of numerous investigations, for reviews see for example 
1,1,0. 



Although the model (|IJ) is relatively simple its classical dynamics is in fact rather 
complicated. In this note we discuss the bosonic sector of ([!]) following the lines of our 
previous paper |T2|. There we have realized that at least at some specials cases, the 



solutions of classical equations of motion 

Xi = [[Xj,Xi],Xj] (2) 

are exponentially unstable, i.e. the system (Q) is stochastic. The appearance of chaos in a 
classical system means that we cannot trust to the ordinary semiclassical analysis of the 
corresponding quantum system. Let us note that a' corrections to the action ([Tp induce 
a stabilization of the classical trajectories Jl3| . 



Here we shall confine ourselves to the simplest version of (]T|) which corresponds to the 
reduction of (2+1) dimensional SU(2) Yang-Mills to (0+1). In the A = gauge we deal 
with eqs.(|2]) for % — 1, 2 and the Gauss law constraint 

[X i ,X i ] = 0. (3) 

High symmetry of the system allows one to reduce the dimension of the phase space. 
Three components of the Gauss law and one more first integral, which we denote by n, 
lead to a four dimensional phase space and the Hamiltonian 

ff= l W+ ^ )+ l W + _^_ + _^_, (4) 

where pf t9 are the momenta conjugated to the f,g. 

In this paper we present analytical and numerical study of the system (|4]). We will 
show that for n large enough the system is integrable and its motion is located in a 
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compact region of configuration space. For n small enough the system exhibits the chaotic 
behavior. Chaotic behavior is a typical feature of systems which one gets as a long-wave 
approximation in a field theory, see for example [|14], [15|, [16], [17 . 



For a better understanding of a role of n-dependent terms in (Mj) we start with the 
toy model governed by the Hamiltonian H with n = and an infinite elastic reflecting 
wall parallel to the g-axis. This model exhibits a chaos-order transition. For the wall far 
enough from the origin the motion is confined to the region p <C / where it admits an 
analytical investigation |TJ|]. We show that this model is integrable in this region. 

As the next approximation (especially valid for g <C / ) we choose the slightly simpli- 
fied version of (f|) with the hyperbolic wall potential: 

We show numerically that it also exhibits chaos-order transition governed by the param- 
eter C, = f^^fr^ an d give some analytical arguments in favor of this result. Note that an 
effect of the extra term 1/2/ 2 to chaotic behavior of a two-dimensional system has been 
discussed in the recent paper fTTjj . 

All this reasoning forces us to conjecture that Hamiltonian (f|) also describes two 
phases depending on the value of n. We compute the Poincare sections for a number 
of characteristic values of n with the energy E and A being fixed. As it expected for 
small n one has the typical stochastic distribution of points and for large n the points are 
distributed along regular lines. 

The paper organized as follows: In Section 2 we present out notation and remind 



the results of [[L2| . Section 3 is devoted to a toy model with an elastic reflecting wall. 
In Section 4 we discuss the model with hyperbolic wall and in Section 5 we present the 
results of numerical calculations. 



2 Notations 



In this section we review the results of |]T2| concerning the appropriate parametrization 
of the configuration space. 

The Lagrangian admits two global continuous symmetries. The SU(2) rotations 

X[ = U + X t U, X' 2 = U + X 2 U, U e SU(2) (2.1) 

yields the conservation of the "angular momentum" 

M = [X 1 ,X 1 ] + [X 2 ,X 2 }. (2.2) 

M = is the Gauss law. The 0(2) subgroup of Lorenz SO(l, 2) 

X l = cos<£>Xi — simp X 2 , 

X 2 = sirup Xi + cos if X 2 , (2.3) 

give one more first integral 

A = Tr(X 1 A 2 -A > 1 X 2 ). (2.4) 
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Hence, it is convenient to parametrize X\ and X 2 as follows: 

X x {t) = V2U + (t) (^f(t) cos9(t) - ^g{t) sin0(*)) U(t), 

X 2 (t) = V2U + (t) (^f(t) sinfl(t) + ^g(t) cos0(t)) U(t), (2.5) 

where f(t), g(t), 6{t) are real functions and U(t) is a SU(2) group element. 

The parametrization ( |2.5| )can be justified as follows. The variables X% and X2 could 
be treated as vectors in the internal isotopic space. At any time they can be rotated to 
belong to some coordinate, say (2,3) plane by using an U(t) 6 SU(2): 

X x = (0, X 2 , X 3 ), X 2 = (0, X 2 , X 3 ), (2.6) 

that fixes U(t) up to rotation around the 1-axis. This rotation could be used to impose 
the following constraint: 

XfXf + X 2 X 3 = 0. (2.7) 
The rotation angle to fulfill (|2.7| ) is 



tan 2% 



{xiy + {x*y - {xiy - {xiy 



The constraint fl2.7|) has a transparent geometrical meaning. Let $ 2 , $ 3 be the fol- 
lowing two- vectors: $ 2 = (X 2 ,Xf) and $ 3 = (Xf,Xf), then (|2~7| ) is the orthogonality 
condition: (<3> 2 ,$ 3 ) = 0. A pair of orthogonal vectors in the plane can be parametrized 
by two radii and one angle (phase), say: 

$ 3 = /(cos 9, sin 9) 

$ 2 = #(-sin#,cos#). (2.8) 

Eqs. ( |2.8|) plus SU{2) rotation give the parametrization ( |2.5| ). Note, that U(l) angular 
momentum N just generates the shifts in 6. 

The main advantage of the coordinate system described above is that four of the six 
Lagrangian equations of motion appear to be nothing but the Noether conservation laws 

M = and N = 0. (2.9) 

Taking into account the Gauss law one gets from Q2.9|): 



2nfg ■ _ n{f 2 + g 2 ) 

ll ~ (/2 _^2)2' 12 ~ U ' /3 - U ' & ~ (/2_^2)2' ^- 10 ^ 



where UU + = I = \Ojlj and n = N. By substituting (|2.10|) into Lagrangian equations for 
/ and g one gets: 

f = n—— T Xfg (2.11) 



(f 2 - 9 2 ' 



n 9 -^Q-Xfg (2,2) 



It is a matter of simple algebra to prove that eqs. fl2.11| ) and ( ]2.12|) are the equations of 
motion following from the Lagrangian 



L = \{f 2 + g 2 )-\f 2 g 2 



n 



n 



%-/) 2 4(g + fr 



(2.13) 



Just this Lagrangian will be the subject of our subsequent analysis. In the particular case 
n = the Lagrangian ( |2.13| ) was the object of intensive study about fifteen years ago in 
the context of long-wave approximation of Yang-Mills theory, this model will be referred 
as the hyperbolyc model. The dynamical system ( |2.13| ) has an additional potential term 
which produce two reflecting walls along the f = g and / = —g axes. The appearance 
of the reflecting walls can crucially change the behavior of the system. To demonstrate 
this in the next section we start with study of the hyperbolic model with elastic reflecting 
wall. 



3 Hyperbolic model with reflecting wall. 



It is well known that the hyperbolic model exhibits a chaotic behavior [|14], [L5| . Equations 
of motion for this system have the form: 



./ 
9 



-Xfg 2 
-\f 2 g 



(3.14) 
(3.15) 



As it is shown in [18fl in the asymptotical regime, where i/ < i one can integrate the 
system of equations ( 13.14] ) and ( |3.15|) by using the Bogolyubov-Krylov method [JTJJ: 



/(*) 

git) 



-a- + /3t + 7 



la t t 

\ A(-Q!y +/3t+j) 6 2 



(3.16) 
(3.17) 



This solution is characterized by four parameters a, [3, 7 and (fo, which are related with 
Couchy's initial data f = /(0), p = /(0) g = g(0) qo = g(0) as: 



fo = 7, Po = P, go 




cosvso, qo 



'2«7 



sin (p 



or more explicitly, 



q 2 o + Vfol 
a = — , cp = arccos 

2 Jo 



foYo 



\ q 2 o + VWo 



Note that the parameter a is the Ehrenfest adiabatic invariant 

f 2 + \f 2 g 2 



a 



2/ 



(3.18) 



(3.19) 



(3.20) 
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In the region / > 0, a is positive and therefore there exists a maximum of coordinate /: 

(3 2 

fmax = 7 + 7T~ (3-21) 

fmax being expressed in terms of dynamical variables is an integral of motion 

f _ (/ 2 + ^ 2 + A/V) rq 

Energy = \(f 2 + g 2 ) + \f 2 g 2 is related to f max as 

E = af max , or / maa; = (3.23) 

Note that a and f max are approximate integrals of motion. 

In the asymptotic region g ~ £/ where one has ^ ~ £v / 2£A and 4 ~ £ 2 A. 

Let us put an elastic reflecting wall at / = I. This means that we consider equation 
( 3.14Q and ( [3.15| ) only for / > I, and g is an arbitrary. We assume that /-component of 



momenta changes a sign upon collision with a wall and g component does not change it. 
For the case of elastic reflecting wall located on / = I there is maximal allowed value of 
g, g max = j\f^jf <^ 1 • The characteristic parameter £ in this case is £ < j2 \f^jf- If one 
consider a trajectory starting from a point on the right of the wall with (<C 1, then the 
trajectory is described rather well by Q3.16 ) and ( 3.17 ). After reflecting the particle moves 



along the trajectory still given be ( |3.16| ) and ( 3.17 ) with new initial data. It is evident 



that the energy and the Ehrenfest invariant are conserved upon a collision of the particle 
with the elastic wall, so the value of maximal deviation is also conserved. Therefore, the 
particle can never reach / = oo. This is basic property of the hyperbolic model with 
reflecting wall located so that the characteristic parameter £ is small enough. In other 
words, we conclude that if we put the reflecting wall so that £ -C 1 then we deal with the 
integrable system and f max is one of its integral of motion. 

4 Model with Hyperbolic Potential. 

In the asymptotic region g ^ f Lagrangian (|2. 13|) has form: 



L =2 if +^-2 f9 2T 2 (4 ' 24) 



as compare with Q2.13D we neglect the o(j2 ) terms. Equations of motion for the Lagrangian 
(|424p are: 



2 n2 



f = -Wf + J5 (4-25) 
g = -Xf 2 g (4.26) 
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We call this model the hyperbolic model with the hyperbolic wall. The equipotential lines 



for ( 4.24j) are defined by: 



A „ o n 



2 



E = 'nf9 2 - ( 4 - 27 ) 

or explicitly: 



2 J y 2/ 2 



1 2EP-n 2 

'-ri-hr- (4 - 28) 



In figure ([!]) and @) we present the form of the potential ( 4.24 ) and ( [2. 13 ), respectively 



and draw corresponding equipotential lines. Simple calculations show that the maximum 
of g is reached at the point / = n/y/E. The minimal accessible value of / is f m i n = 
n/\ / 2E. So, we deal with the potential for which equipotential lines go to infinity along 
the /-axis. Maximal value of g/f in terms of E, n and A is: 

^(fL^yivx (4 - 29) 

The condition £ -C 1 guaranties that the system is always in the asymptotic region g ~ £/ 
with £ < 1. 



We integrate this system asymptotically by using Bogolyubov-Krylov method |19 
According to this method one has to take 



9{t) = yXf(tj cosipit) (4 - 30) 



with some constant a, then integrate the equation 



dt 2 p 
Equation ( f4.31|) can be easily integrated, 



d 2 f n 2 , . 

+ a ~ « = ( 4 - 31 ) 



'2a(f max f) , 2 f2 (fmax f 2 

J J max 

= y/\ J f(t)dt + <p (4.33) 

Here constant of integration is chosen so that at t — the /-coordinate takes its maximal 
value. Energy E and maximal deviation f max are related as follows: 

n 2 

E = af max + -— (4.34) 

J max 

One can invert this algebraic equation and get f m ax as function of dynamical variables of 
our system fl4.24p . 



As result we suggest that in close analogy with previous model for ( <^ 1 the sys- 
tem becomes integrable and f max is integral of motion. We shall justify this conjecture 
numerically. 
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5 Numeric calculations. 



In this section we investigate the existence of chaos-order transition for the hyperbolic 
model with reflecting wall, model with hyperbolic wall ( [4.24 ) and system fl2.13| ) numeri- 
cally 

At first we test the Poincare sections for the model with a reflecting wall. The con- 
servation of energy restricts any trajectory of the four- dimensional phase space to a 
three-dimensional energy shell. At a given energy any additional constraint defines a 
two-dimensional surface in the phase space, which is called the Poincare section It is con- 
venient to take a constraint g = 0. All crossections of a trajectory with the surface are 
marked by points on the (/, pj)-plane. On each figure we plot a Poincare section for a set 
of trajectories to show that behavior of the system does not depend on the initial data. 
Any trajectory was integrated as long as the program guarantees that the deviation of 
the energy is less then 0.1%. Different colors correspond to different trajectories with 
fixed parameters and random initial data. Chaotic motion is characterized by a set of 
randomly distributed points. Regular trajectories are depicted by dotted curves. 

In Figures (|3]), (|) and (|5|) we plot Poincare sections for the different values of the 
reflecting wall coordinate / and with the same energy E = 1. The pictures show that for 
small values of the parameter / chaotic region is located near the wall. For large values 
of the parameter I the points arrange into closed dotted curves. 

An important characteristic of a dynamical system is the Lyapunov exponent f](t). It 
has a positive limit: (lim^oo r] > 0) for a chaotic system and zero limit: (lim^oo r] = 0) 
for a regular one. The calculations for the model with the hyperbolic wall ( 4.24 ) were 
performed for different values of ( (|4.29| ). By changing the energy E with A and n being 
fixed. We vary the parameter (. The program starts with random initial data with 
given energy. The program calculates the coordinates / and g, the energy, f max and the 
Lyapunov exponent. Typical results for ( = 1 and ( <C 1 are shown on Figures (||) and 
(0) respectively. One can see that for ( = 1 (Fig. || white curve) Lyapunov exponent has 
a positive limit. For ( <C 1 the Lyapunov exponent goes to zero (Fig. [7|, white curve). 
Parameter f max (blue curve) for small ( does not change with the time. Energy (red 
line) shows that program works perfectly and energy does conserve through the whole 
calculation time. The numerical calculations show that for ( = 1 the system ( |4.24j ) is 
stochastic and for £ <C 1 it is the integrable one, that confirm the analytical results of 
sect. 4. 

Now we turn to the main model ( [2.13| ). The dynamics of ( |2.13| ) was analyzed in |L2| 
for relatively small (in units of E) values of n and we have found it to be fully chaotic. 
For the relatively large values of n the motion is confined into the region ( <C 1 (see Fig. 
D. In this region the model with hyperbolic wall Q4.24|) gives a good approximation to 
the main model ( p.l3j ). We conjecture the main model also to be integrable for the large 
n. In favor of this conjecture in Figures (|8]), (|9]) and (|K]) we plot the Poincare sections 
for ( |2.13D with different values of parameter n and fixed energy E = 1. For large values of 
parameter n there are only regular closed orbits. We also test the Lyapunov exponent and 
obtained the pictures quite similar to the ones for the model (|4.24|). Therefore we see that 



the system exhibit the chaos-order transition governed by the characteristic parameter 
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Figure 4: Reflecting wall potential, I = 10, E = 1 



11 



re 5: Reflecting wall potential, I = 50, E 




Figure 6: Chaotic phase 
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Figure 7: Regular phase 




Figure 8: Full potential, n — 1, e = 1 
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